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ABSTRACT 

Although variations in elemental abundance ratios in the Milky Way Galaxy certainly exist, 
details remain uncertain, particularly in the inner Galaxy, where stars and H II regions in the 
Galactic plane arc obscured optically. In this paper we revisit two previously studied, inner 
Galaxy H II regions: G333.6— 0.2 and W43. We observed three new positions in G333.6— 0.2 
with the Kuiper Airborne Observatory and reobserved the central position with the Infrared 
Space Observatory's Long Wavelength Spectrometer in far-infrared lines of S ++ , N ++ , N + , and 
++ . We also added the N + lines at 122 and 205 to the suite of lines measured in W43 by 
Simpson et al. (1995). The measured electron densities range from ~ 40 to over 4000 cm -3 in 
a single H II region, indicating that abundance analyses must consider density variations, since 
the critical densities of the observed lines range from 40 to 9000 cm -3 . We propose a method to 
handle density variations and make new estimates of the S/H and N/H abundance ratios. We 
find that our sulfur abundance estimates for G333.6— 0.2 and W43 agree with the S/H abundance 
ratios expected for the S/H abundance gradient previously reported by Simpson et al., with the 
S/H values revised to be smaller owing to changes in collisional excitation cross sections. The 
estimated N/H, S/H, and N/S ratios are the most reliable because of their small corrections for 
unseen ionization states (< 10%). The estimated N/S ratios for the two sources are smaller than 
what would be calculated from the N/H and S/H ratios in our previous paper. We compute 
models of the two H II regions to estimate corrections for the other unseen ionization states. We 
find, with large uncertainties, that oxygen does not have a high abundance, with the result that 
the N/O ratio is as high (~ 0.35) as previously reported. The reasons for the uncertainty in 
the ionization corrections for oxygen are both the non-uniqueness of the H II region models and 
the sensitivity of these models to different input atomic data and stellar atmosphere models. We 
discuss these predictions and conclude that only a few of the latest models adequately reproduce 
H II region observations, including the well-known, relatively-large observed Ne ++ /0 ++ ratios 
in low- and moderate-excitation H II regions. 
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1. Introduction 

The presence of radial abundance gradients in 
the plane of the Milky Way is now an established 
fact and is seen in both stars and gaseous nebu- 
lae (see, e.g., Henry & Worthey 1999; Rolleston 
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et al. 2000; and references in both papers). The 
cause is generally thought to be due to the for- 
mation history and subsequent evolution of the 
Galaxy; thus the observed gradients are a ma- 
jor tool for understanding this history (see, e.g., 
Hou, Prantzos, & Boissier 2000; Chiappini, Mat- 
teucci, & Romano 2001; Chiappini, Romano, & 
Matteucci 2003, and references therein). The ba- 
sic idea is that the inner Galaxy formed before the 
outer Galaxy and the higher molecular gas den- 
sity in the inner Galaxy produces a higher star 
formation rate. The result is a greater return to 
the interstellar medium (ISM) in the inner Galaxy 
of both "primary" alpha elements from massive 
star supernovae and "secondary" elements like ni- 
trogen. Secondary nitrogen is produced by CNO 
burning of already existing carbon and oxygen 
in intermediate-mass stars and subsequently re- 
turned to the ISM through mass loss. However, 
both the wide range of uncertain input param- 
eters to chemical evolution models (Pagel 2001) 
and the uncertain details of the abundance varia- 
tion of each element, primary and secondary, con- 
tribute to our less-than-complcte understanding of 
the formation and evolution of the Milky Way. 
In this paper we try to improve our knowledge 
of the abundances in inner Galaxy H II regions 
and also test the methods that have been used to 
date to determine the abundance ratios from low- 
excitation, optically obscured H II regions that are 
needed for abundance gradient calculations. 

H II regions in the inner Galaxy (galactocen- 
tric radius Rq < 6 kpc, where R Q — 8 kpc) are 
not accessible to optical observers because of inter- 
stellar extinction. Consequently, the abundances 
must be determined from far-infrared (FIR) lines, 
such as those of N+ at 122 and 205 /im, N++ at 
57 /im, 0++ at 52 and 88 /im, Nc+ at 12.8 /im, 
Ne++ at 15.5 and 36 /im, S++ at 18.7 and 33.5 
/im, and S 3+ at 10.5 /im. Large surveys include 
those of Simpson et al. (1995a, hereafter SCREH) 
and Afflerbach ct al. (1997), who observed the 
FIR lines (wavelengths > 17 /im) with the Kuiper 
Airborne Observatory (KAO) in 18 inner Galaxy 
H II regions, and Martin-Hernandez et al. (2002a; 
2003), who observed 13 inner Galaxy H II regions 
with the Infrared Space Observatory (ISO). The 
KAO observations included FIR measurements of 
the lines with wavelength > 18 /im, whereas the 
ISO observations included the mid-infrared lines, 



measured with ISO's Short Wavelength Spectrom- 
eter, SWS, but not the [N II] 205 /im line. 

Two of the abundance gradients that are most 
needed for galactic chemical evolution studies are 
the O/H and N/H ratios, representing primary 
and secondary elements, respectively. In optically 
obscured H II regions the only ionization state of 
oxygen available for study is ++ , which means 
that corrections are required for the unseen + 
ions. The situation is better for nitrogen, as there 
are [N II] lines as well as [N III] lines in the FIR 
wavelength range. Unfortunately, the [N II] 122 
/im line is difficult to measure for airborne observa- 
tions because of atmospheric absorption, and the 
[N II] 205 /im line is difficult because of poor detec- 
tor sensitivity. In addition, the [N II] lines are dif- 
ficult to interpret because they are collisionally de- 
excited at much lower densities than the lines from 
the other ions. Consequently, the papers to date 
have used only the [N III] and [O III] lines to esti- 
mate N/O and the secondary /primary abundance 
ratio in optically obscured H II regions (SCREH; 
Afflerbach et al. 1997; Martin-Hernandez et al. 
2002a). This is particularly troubling since a ma- 
jority of the nitrogen and an even larger frac- 
tion of the oxygen are singly ionized in the high- 
metallicity H II regions of the inner Galaxy (note 
that in the outer Galaxy, the lower metallicity H II 
regions can be much more highly excited and thus 
the N ++ /0 ++ ratio is often a good approximation 
to N/O, Rudolph et al. 1997). On the other hand, 
a majority of the primary element sulfur is doubly 
ionized in H II regions. Here we devise a method of 
analyzing the density-sensitive [N II] lines to more 
accurately estimate nitrogen abundances with re- 
spect to sulfur without the need for large correc- 
tions for unseen ionization states. Using the ob- 
served N++/N+ ratio and calculated H II region 
models, we then estimate the ionization correction 
factors (ICFs) for oxygen. 

Two of the most luminous H II regions in the 
Galaxy are G333.6— 0.2 at Rq ~ 5.5 kpc (Wilson 
et al. 1970) or 5.6 kpc (Caswell & Haynes 1987) 
and W43 (G30.8-0.0) at R G ~ 4.2 kpc (Reifen- 
stcin ct al. 1970) for R Q = 8.0 kpc. The radio 
fluxes from these objects require ionizing stellar 
luminosities of 10 50 and 2.9 x 10 50 photons s _1 , 
respectively (see Rubin 1968b). Both H II re- 
gions are located in the Scutum-Crux spiral arm 
(see Vallee 2002) in the Milky Way molecular ring. 
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Both are excited by clusters of massive stars, ob- 
served at near-infrared (NIR) wavelengths. The 
cluster exciting G333.6— 0.2 contains stars that are 
very young, embedded in dust (Blum et al. 2002), 
in keeping with the youthful appearance of this 
extremely dense, compact H II region. On the 
other hand, the cluster exciting W43 must be older 
and more evolved because it contains a Wolf-Rayet 
star in addition to O giants or supcrgiants (Blum, 
Damineli, & Conti 1999; Cotera & Simpson 1997), 
although there are still numerous stars forming in 
the associated molecular clouds (Motte, Schilke, 
& Lis 2003). This is also consistent with the more 
mature appearance of W43, where the ionized gas 
is well separated from the exciting star cluster. 

In this study, we report measurements of the 
[N II] lines in the three positions observed in W43 
by SCREH and numerous lines in three additional 
positions in G333.6— 0.2 north of the center po- 
sition observed by Colgan et al. (1993) with the 
KAO; we also report measurements of the central 
position in G333.6— 0.2 with ISO's Long Wave- 
length Spectrometer (LWS). In Section 2 we de- 
scribe the observations, in Section 3 we discuss 
analysis techniques that are used to determine 
abundances in both low- and high-density H II re- 
gions, in Section 4 we model the two H II regions 
to determine the ICFs needed for those elements 
without measured abundances of the singly ion- 
ized atoms, and in Section 5 we discuss models in 
general and give our conclusions. By measuring 
both N+ and N ++ , we hope to avoid the need for 
large (and uncertain) ICFs for nitrogen and thus 
improve the reliability of the secondary /primary 
abundance ratio measurement in the Milky Way 
molecular-ring H II regions. 

2. Observations 

Both G333.6— 0.2 and W43 are extended com- 
pared to our 45-60" KAO beam and 80" LWS 
beam, which means that integrating over a con- 
stant density model to calculate ICFs is not ade- 
quate. Instead, we measured the fluxes at a num- 
ber of positions at varying distances from the ex- 
citing stars to compare with the radial distribu- 
tions of line flux computed from models to deter- 
mine the source excitation and abundances. 



2.1. KAO Observations 

The observations were made with the Cryogenic 
Grating Spectrometer (CGS, Ericksonet al. 1985) 
on NASA's 0.91 m KAO telescope. Data acqui- 
sition and reduction techniques are described by 
Colgan et al. (1993) and SCREH. For the present 
observations, Saturn was used for calibrating the 
absolute fluxes and determining the relative detec- 
tor responses. The observed line and continuum 
fluxes for all positions are given in Table 1. 

G333.6— 0.2 was measured at its center and in 
three offset positions. The data from the center 
position, C (RA 16 h 22 m 9!3, Dec -50°6'5", both 
J2000), were taken with a 45" beam and were pre- 
viously reported by Colgan et al. (1993). The 
three North positions (Nl, N2, and N3 at 1', 2', 
and 3' north of G333.6— 0.2C) were measured on 
1993 May 11, 13, and 20, flying from Christchurch, 
New Zealand. For these positions, a 60" beam was 
used for all lines to provide good ratios to the [N II] 
205 /im line, at whose wavelength the KAO has a 
diffraction size of 1.22X/D = 57". The chopper 
angle was approximately east-west (position angle 
~ 90°) and the throw was ~ 8.5'. Because the 
[N II] 122 /im line is on the edge of a strong tel- 
luric H2O line, the observations of positions Nl, 
N2, and the [N II] 205 line in Position N3 
were taken at or near 43,000 feet altitude, instead 
of the more usual 39,000-41,000 feet, where the 
other lines in Position N3 were measured. Cor- 
rections for telluric H 2 absorption were made 
for both source and calibrator — we estimate the 
zenith H2O overburden to range from 4.4 to 5.2 
/jm at 43,000 feet on May 11 and 13, and 3.9 ^m 
at 41,000 feet on May 20. 

The three positions of W43 that had pre- 
viously been studied by SCREH (W43N at 
RA 18 h 47 m 36!8, Dec -1°55'29"; W43C at RA 
18 h 47 m 36!0, Dec -1°56'39"; W43S at RA 18 h 47 m 38?5, 
Dec -1°57'34", all J2000) in the [N III], [O III], 
and [S III] lines were measured in the [N II] lines 
on 1992 July 1, flying from NASA Ames. The 122 
jim line was measured in both 45" and 60" beams 
and the 205 /im line was measured in a 60" beam. 
The chopper angle was approximately cast-west, 
the throw was ~ 6', and the zenith telluric H2O 
overburden was ~ 5 (im at 43,000 feet and ~ 2 
/j,m at 45,000 feet. We estimate that there are 
extra measurement uncertainties of order 20% in 



3 



addition to the statistical errors given in Table 1 
for the W43 [N II] 122 lines. These could 
result from dumpiness, pointing differences, or 
the effects of poor correction for deep atmospheric 
absorption. The analyses of the various measure- 
ments produce different abundance measures, as 
will be tabulated in the next sections. 

2.2. ISO Observations 

G333.6— 0.2 was measured with the ISO Long 
Wavelength Spectrometer in LWS01 mode on 1997 
March 20 (TDT No. 49001002, PI R. Rubin), ex- 
posure time 1054 s. The line strengths for the 
[O III] 52 and 88 ^m and [N III] 57 fim lines are 
given in Table 1 (the [O I] 63 ^im line is also strong, 
but the [N II] 122 ^m line was not detected owing 
to ISO's low spectral resolution in LWS01 mode 
with the LW detectors). 

2.3. Extinction Corrections 

We corrected the measured line fluxes for ex- 
tinction using the infrared extinction function of 
Li & Draine (2001). This function has a very simi- 
lar wavelength dependence to the extinction func- 
tion used by SCREH when normalized to Li & 
Drainc's extinction maximum at ~ 9.5 /im. Thus 
the FIR optical depths are characterized as being 
proportional to the maximum optical depth in the 
9.5 to 10 /um region, T9.7. Based on new infor- 
mation in the literature, we have updated these 
values from SCREH, with W43 having more ex- 
tinction and G333.6— 0.2 having less. 

For three stars in the central cluster of W43, 
Blum et al. (1999) estimated that Ak = 3.55, 
3.52, and 3.63. Using the extinction law of Li 
& Draine (2001), we find that Ak = 3.55 cor- 
responds to T9.7 = 2.74. Thus we increase the val- 
ues of T9.7 used by SCREH by 37% to t 9 . 7 = 2.74, 
2.74, and 3.63 for W43N, W43C, and W43S, re- 
spectively. 

For G333.6-0.2, Fujiyoshi ct al. (2001) inferred 
an average value of T9.7 = 1.5 from their mid-IR 
spectropolarization measurements. Consequently, 
we use 79.7 = L5 for the center position instead 
of SCREH's T9.7 = 2.0. (This value corresponds 
to An a = 14 mag; a lower extinction yet for the 
center of G333.6— 0.2 might be predicted, since the 
nebula was detected in Ha by Churms et al. 1974.) 
We use T9.7 = 1.5 for the North positions as well, 



since there are no significant variations in the in- 
tegrated H I optical depths measured by Forster 
et al. (1987) at the C, Nl, and N2 positions. 

Following SCREH, an extinction uncertainty 
estimate of 20% (la) is propagated through all 
subsequent data analysis errors. 

3. Analytical Approach 

Because recombination, brcmsstrahlung, and 
collisional excitation rates are all proportional to 
the electron density squared, iV|, it has long 
been the practice to treat all abundance analyses 
in H II regions and planetary nebulae (PNe) as 
though the surface brightness is, in fact, propor- 
tional to the emission measure, EM = J N^dl, 
where I is the line-of-sight path through the neb- 
ula. This is valid so long as N e is substantially 
less than the critical density, N cr n, defined as 
the density at which the upper level of a tran- 
sition is depopulated equally by radiative transi- 
tions and by collisions (e.g., Osterbrock 1989; Ru- 
bin 1989). Corrections for collisional de-excitation 
can be made if N e is known and approximately 
constant. However, constant densities arc un- 
likely to occur in real nebulae. The effect of den- 
sity variations is that the kernels in the integral 
over the nebular volume have different weighting 
functions for ions with different amounts of colli- 
sional de-excitation. Thus the flux ratios of such 
ions depend on these variable weighting functions 
(see, e.g., Rubin 1989). Here we consider the 
calculation of ratios of column densities of ions 
whose emission is proportional to N e , as well as to 
Af?. The former is appropriate when N e 3> N crit . 
These two methods give identical results for ion 
ratios when the density is constant over the whole 
observed volume. We will demonstrate the effects 
of density variations via a simple example and use 
the results to estimate how to handle the transi- 
tion case where N e ~ N cr it- 

3.1. Low Density H II Regions 

For a given ion, a "low density" H II region 
is one where N e <C N crit and a "high density" 
H II region is one where N e ^> N cr it\ thus any 
particular H II region can have both low density 
and high density volumes, depending on the value 
of N cr i t for the transition being studied. 

In the optically thin case, the observed line flux 
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is F — J jdldfl where j is the volume emissiv- 
ity divided by 4-zr and dldQ is the volume ele- 
ment of the H II region in the telescope beam, 
f2. We define the normalized volume cmissivity 
e = 4nj/NiN e (Rubin et al. 1994; SCREH) where 
iVj is the density of atoms in the particular ion- 
ization state and the average 7V e -weighted ionic 
density as < N t >= / NiN e dldQ/ J N e dldQ. For 
the low density regime where N e -C N crit , e has 
very little dependence on N e , and the line surface 
brightnesses (or fluxes) are proportional to J N^dl 
(or jN^dldQ). 

Thus for this case we can write 

F = [ -^-NiN e dldn w 4~ < > [ N e dldtt 
J 4?r 4tt J 

(1) 

and the ionic abundance ratio is 

<N il > = Fi/ei . , 

<N l2 > F 2 /e 2 - [ > 

As N e approaches N cr n, collisional de-excitation 
of the line's upper level becomes more and more 
important, and the ratio of two lines from the same 
ionic species iVj with different N cr i t can be used 
to estimate the electron density (e.g., Rubin 1989; 
Fig. la). 

3.2. High Density H II Regions 

At high densities, the fractional population of 
the upper level approaches that of local thermo- 
dynamic equilibrium (LTE) and is a function of 
electron temperature, T e , only. In this regime, the 
line surface brightness is proportional to the col- 
umn density, CD = J Nidi, the number of ions in 
the line of sight and 

F = J ^-Nididn « J N^i, (3) 

where the quantity N e e is the emissivity per ion 
and has very little density dependence so long as 
there are no low density regions along the line of 
sight (no large density variations). Then for iden- 
tical beam areas f2, the column density ratio is 

CDi = F 1 /(N ei e 1 ) 

CD 2 F 2 /(N e2 e 2 )- [ ' 

In fact, for abundance ratio calculations, one really 
wants to measure the actual number of ions in the 



line of sight, and not the number of ions weighted 
by some other, strongly varying parameter such as 
the electron density (as in equation (4)). Thus it is 
desirable to estimate column densities for gaseous 
nebulae whenever feasible. 

3.3. Density Variations: An Illustrative 
Model 

The presence of density variations has a number 
of effects on the abundances that one infers from 
FIR line measurements. We demonstrate these ef- 
fects by the use of simple slab models, with an ex- 
ponentially decreasing density N e = N e~ x from a 
location, xq, in the interior of the slab to an outer 
location, x max . This density function was chosen 
because it gives similar weighting to both high and 
low density regimes (other types of models, e.g., 
Rubin 1989, could be used to illustrate the same 
effects). For simplicity, the example assumes no 
ionization stratification, N{ = N e , and T e = 7000 
K; thus the only effects on the calculated volume 
emissivities are those of the variation in density. 
For all models, a very large value of x max was 
chosen such that N e — 0.01 cm~ 3 . Each model 
had a different value of xo and A*o, ranging from 
No = 1 cm~ 3 to A*o = 10 6 cm~ 3 . The intensities, 
/, emitted in each line were calculated from 

/ = ^-NiN e dx. (5) 

Figure la shows how the electron density, N est , 
can be estimated using calculated volume emis- 
sivity ratios from the observed [S III] 19/33 /xm, 
[O III] 52/88 yum, and [N II] 122/205 ^m line ra- 
tios, with the [N II] ratio being sensitive to low 
densities and the [S III] ratio sensitive only to 
much higher densities. Values of N est were com- 
puted from the ratios of the line intensities pro- 
duced by the models; ratios of these estimated 
densities to the maximum density of each model, 
TVo, are plotted in Figure lb. This figure shows 
that in the presence of density variations, even 
when the ions coincide exactly in space, one does 
not infer the same density from the three ratios: 
N est (N II) < N est (0 III) < A est (S III). 

To understand the effects of density variations 
on estimated abundances, we calculate the esti- 
mated numbers of ions in the volume from the 
model line intensities and N est and compare these 
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estimates with the true numbers of ions found 
by integrating N e or N% over the model line of 
sight. These comparison ratios we call "abundance 
measures" . Numerically, we define the abundance 
measure for the low-density case where the line in- 
tensities are proportional to the emission measure 
as 

I 



AM EM = 



NN e dx 



(6) 



and the abundance measure for the high-density 
case where the line intensities are proportional to 
the column density as 



AM C d = 



(N est e) 



Ndx 



(7) 



AM em and AMqd are plotted in Figure 2 for 
each ion as functions of N est . Both abundance 
measures contain fewer ions than the input num- 
ber, but the regimes of reasonable validity for the 
two measures are different and depend on N crit . 
Fortunately, much of the discrepancy with the N e 
or TVg volume integrals cancels when we take ratios 
of the EM or CD abundance measures (equivalent 
to equations (2) and (4), respectively), as shown 
in Figure 3. These plots suggest an approach for 
the intermediate cases where N es t ~ N cr it and 
there are measurements of N e from (at least) two 
different density-sensitive line pairs: N ei and N e2 . 

Using Fig. 3, we define three regimes as fol- 
lows: The low-density regime occurs when the es- 
timated density is lower than the lowest N cr it for 
the lines involved (see Fig. 2 for N cr i t ) — here 
we use the equations for emissivity oc since 
the ionic ratios weighted by are close to the 
model input ratio (unity for the plots). The high- 
density regime occurs where the column density 
ratio is closer to the model input ratio than the 
ionic ratios weighted by N% — this occurs when 
the estimated density is higher than the larger 
N cr u of the more density-sensitive line pair. At 
intermediate densities, both the low-density and 
high-density equations give abundances discrepant 
from the model input. Considering this, in the in- 
termediate regime we average the two cases and 
quote an uncertainty equal to the greater of one 
half the difference between the two abundance ra- 
tios and the statistical error. For the low- and 
high-density regimes, we will use the equations for 



the ionic ratios weighted by or the column 
density ratios, respectively. Unfortunately, it is 
extremely difficult to get a good estimate of the 
column density of H + because the emissivities of 
H lines in H II regions and PNe are always pro- 
portional to iVg and the properly weighted den- 
sity is impossible to determine. Because of this, 
we are not able to measure abundance ratios with 
respect to hydrogen for all the ions; instead we 
will obtain standard "ionic ratios" with respect to 
H+ only for those ions whose lines have high N crit 
compared to the densities estimated from the ob- 
served line pairs. We will calculate both ionic ra- 
tios and column density ratios for the rest of the 
heavy element ions observed at FIR wavelengths 
and consider the density regime according to the 
line's N crit for estimating the final abundance ra- 
tio. 

A significant difference between this method 
and previous abundance determinations from FIR 
lines is the increase in the size of the errors that 
occur when more than one density-sensitive line 
pair is measured in a high-density, compact H II 
region. The smaller statistical errors that one ob- 
tains when only one density line pair is measured 
are an illusion. It is clear that measuring multiple 
density-sensitive line pairs is not just desirable but 
is, in fact, required to obtain reliable abundance 
measurements from FIR lines (even though they 
may have large assigned errors). We note that a 
similar situation occurs in optical abundance de- 
terminations from forbidden lines in that T e must 
be inferred from various line ratios to obtain re- 
liable abundances (e.g., Kennicutt, Bresolin, & 
Garnctt 2003). 

4. G333.6 0.2 and W43 

4.1. Analytical Abundance Estimates 

The complete set of derived N e values and 
abundance ratios is given in Table 2. The val- 
ues of T e used were 6200 K (Caswell & Haynes 
1987) and 6500 K (Wink, Wilson, & Bieging 1982) 
for G333.6— 0.2 and W43, respectively (mid-and 
FIR lines are not sensitive to uncertainties in T e ). 
The radio fluxes for G333.6— 0.2 were estimated 
by integrating over the model described in the 
next section at the appropriate distances from the 
center (Figure 4) and the radio fluxes for W43 
were taken from Lester et al. (1985) (the NRAO 
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VLA Sky Survey, Condon et al. 1998, has lower 
fluxes, probably because the bright clumps in W43 
are somewhat optically thick at 1.4 GHz, and the 
"snapshot" observation of Balser et al. 1995 is 
missing extended flux). The differences between 
the ionic abundance results presented here and 
in SCREH for the KAO data are due to differ- 
ent extinction corrections (§2.3), different radio 
fluxes for G333.6— 0.2, and different collisional ex- 
citation cross sections for [S III] (Tayal & Gupta 
1999). The differences between the ISO and KAO 
abundances with respect to H + probably indicate 
the uncertainties in the radio flux estimates. The 
N++/0++ ratios are the same to within the errors 
(3cr), although some of the difference may be due 
to different beam sizes or pointing errors. 

Estimated abundance ratios for N/S and S/H 
are given in Table 3, where we assumed that 
the ionization fraction for sulfur <S++/S>= 0.9 
(a typical value for the low-to-moderate excita- 
tion and high-metallicity models in Rubin 1985, 
SCREH, and the models computed in §5) and that 
the ionization fraction of <N 3+ /N> is so small 
that N 3+ can be neglected. The N/H ratio was de- 
rived by multiplying N/S times S/H. Abundance 
estimates for the other elements are more uncer- 
tain because they require the use of models to es- 
timate the necessary ICFs. 

4.2. Models of G333.6 0.2 and W43 

We computed H II region models using the 
program NEBULA (Rubin 1968a; Rubin 1985; 
SCREH). We have upgraded NEBULA since the 
discussion in SCREH in the following ways: 

(1) The photoionization and recombination 
cross sections were extensively updated to include 
results from the Opacity Project (Seaton et al. 
1992). Fits to photoionization cross sections by 
Verner et al. (1996) were used for photoionization 
of C+ C++, N°, N+, N++, 0°, 0+, Ne°, Nc+, S+, 
S++, S 3 +, and Si+. The changes to the coefficients 
for ionization and recombination of iron are dis- 
cussed by Rodriguez & Rubin (2004). (Changes 
to other elements were omitted because of the 
lack of corresponding Opacity Project recombina- 
tion coefficients.) Total recombination coefficients 
(including diclcctronic recombination) were taken 
from the tabulations of Nahar & Pradhan (1997) 
for C+, C++, and N++; Kisielius & Storey (2002) 
for N+; Nahar (1999) for O and 0+; Kisielius et 



al. (1998) for Ne+; Nahar (2000) for S++; and Na- 
har (1995) for S+ and Si+. Charge exchange rate 
coefficients for reactions with H+ were taken from 
the tabulation of Kingdon & Ferland (1996) with 
updates from the Oak Ridge National Laboratory 
web site "Charge Transfer Database" 1 and coeffi- 
cients for reactions with He+ were taken from J. 
G. Wang ct al. (2004, in preparation), with the 
data available from the same Oak Ridge web site. 

(2) The collisional excitation cross sections now 
include N+ and 0++ from Lcnnon & Burke (1999), 
N++ from Blum & Pradhan (1992), Ne+ from 
Saraph & Tully (1994), Nc++ from McLaughlin 
& Bell (2000), S++ from Tayal & Gupta (1999), 
S 3 + from Saraph & Storey (1999) or Tayal (2000) 
for T e > 10, 000 K, Ar+ from Pclcn & Berring- 
ton (1995), and Ar++ from Galavfs, Mendoza, & 
Zeippen (1995, 1998). 

(3) The diffuse radiation field calculation (Ru- 
bin 1968a) was modified to include calculations 
at all radial grid positions instead of interpolation 
between selected positions. Up to 20 iterations 
were used if needed to obtain uniform convergence 
(most positions converged within 10 iterations). 

(4) Non-LTE stellar atmospheres with winds 
from Pauldrach, Hoffmann, & Lennon (2001) were 
used. Smith, Norris, & Crowther (2002) have also 
computed model atmospheres using Pauldrach's 
code; however, we did not use Smith ct al.'s mod- 
els because of the coarseness of their wavelength 
grid. 

We used the set of model atmospheres, which 
were downloaded from Pauldrach's web site 2 , con- 
sisting of models for Solar composition, dwarfs (D) 
and supergiants (S), and T e s = 30, 35, 40, and 45 
kK. Although G333.6-0.2 and W43 have metal- 
licitics somewhat larger than that of the Orion 
Nebula, the Orion Nebula abundances (e.g., Simp- 
son et al. 1998) are sufficiently below the "Solar" 
composition used in the stellar atmosphere mod- 
els that G333.6— 0.2 and W43 are probably closer 
in metallicity to the "Solar" composition models 
than to anything greatly different. Stellar atmo- 
spheres with higher metallicities have more winds, 
which produce more extreme ultraviolet (EUV) 
flux, but also more line blanketing (less EUV); the 
difference is probably mainly for frequencies > 4 

1 http://www-cfadc.phy.ornl.gov/astro/ps/data/ 
2 http://www.usm. uni-muenchen.de/pcoplc/adi/adi. html 
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Ryd, which are not important for H II regions. 
We computed models for W43 using both D and S 
models from Pauldrach et al. (2001) and also the 
non-LTE, plane-parallel model atmospheres with- 
out winds calculated by Lanz & Hubeny (2003) 3 . 
We used only the dwarf star model atmospheres 
for G333.6— 0.2 because it is clearly excited by a 
very young star cluster (Blum et al. 2001). 

One result of the changes to the nebular atomic 
physics data is that the volume-averaged ion- 
ization fractions (SCREH; Rubin et al. 1994) 
for <0++/0> and <Nc++/Ne> predicted by 
the models are larger than previous predic- 
tions, whereas the fractions for <S ++ /S> and 
<N ++ /N> are approximately the same. The dif- 
ferences are similar to typical uncertainties in ion- 
ization ratios computed from observational data 
(< 30%), which would make testing difficult un- 
less the data sample was large and of very high 
quality. However, the choice of stellar atmosphere 
model makes a much larger difference in the ion- 
ization fractions and ratios thereof; this choice 
has serious implications for estimates of T c s and 
ICFs. We will compare various volume-averaged 
ionization fractions computed for different model 
H II regions as a function of stellar atmosphere 
model and discuss the need for winds further in 
§5.1. 

4.2.1. G333. 6-0.2 

The initial model is that of a spherically sym- 
metric, core-envelope H II region (Rubin, Hollcn- 
bach, & Erickson 1993), with density decreasing 
from high values in the core to low values in the 
outer envelope, as suggested by the estimated den- 
sities of Table 2. Our initial density distribution 
for the G333.6— 0.2 model computed here was de- 
rived by inverting the radio brightness tempera- 
ture distribution, since the radio surface bright- 
ness is proportional to J N^dl. The radio bright- 
ness temperature distribution of G333.6— 0.2, of 
necessity, has been estimated by combining a va- 
riety of different measurements, none of which in- 
cludes the whole H II region at a frequency high 
enough to insure that the emission is everywhere 
optically thin. The highest spatial resolution ra- 
dio maps in the literature consist of the following 
maps: (1) a Parkes 64-m telescope map at 8.87 

3 Available from http://tlusty.gsfc.nasa.gov 



GHz with half-power beamwidth, HPBW = 2.5', 
made by McGee, Newton, & Butler (1979); be- 
cause the map was made with a single dish, it 
covers the full spatial extent of G333.6— 0.2 to a 
radius of > 15'; (2) a Fleurs Synthesis Telescope 
map at 1.415 GHz with HPBW - 50" by Retal- 
lack & Goss (1980); only the inner 5' could be 
mapped, and the core is optically thick; (3) an 
Australia Telescope Compact Array map at 3.4 
cm (8.8 GHz) with resolution 1.6" by Fujiyoshi 
(1999); only the inner 13" could be mapped, and 
the peak brightness temperature of ~ 4400 K in- 
dicates that the central 2" — 3" is not completely 
optically thin (T e ~ 6200 K, Caswell & Haynes 
1987). However, this last map is extremely im- 
portant because it shows just how compact and 
bright the central core of G333.6— 0.2 really is. A 
compact bright core is also required by the NIR 
maps and the density-sensitive FIR line ratios. 

Still lacking is a measurement of the transition 
region between the inner core and the envelope. 
There is, however, a map of the warm dust emis- 
sion made by the Midcourse Space Experiment 
(MSX; Price et al. 2001) at 21 ^m with 19" reso- 
lution. 4 Although there is no certainty that warm 
dust has the same spatial structure as ionized gas 
(and it certainly does not have constant tempera- 
ture), in the absence of an alternative we used it 
to interpolate between the core and envelope radio 
maps. 

The total radio flux was taken from the mea- 
surements at the highest observed frequencies, v, 
because G333.6— 0.2 shows evidence of being at 
least partially optically thick at frequencies < 10 
GHz. We used the integrated flux of 93.8 Jy at 
14.7 GHz measured by McGee & Newton (1981). 
The final estimated radio brightness distribution 
was renormalized to produce this total integrated 
radio flux. The estimated radio brightness tem- 
perature distribution, scaled oc u~ 21 as though it 
were optically thin at 5 GHz, and derived electron 
densities are shown in Figure 4. The total radio 
flux requires a stellar ionizing photon luminosity 
of 10 50 photons s _1 . There may also be, of course, 
additional EUV photons that do not contribute to 
the H II region photoionization because they are 
absorbed by dust. Table 4 lists the relevant IR 
and radio observations, corrected for extinction, 

4 Available at http://irsa.ipac.caltech.edu 
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including those mid-IR observations from the lit- 
erature which include most of the core. 

In addition to matching the radio brightness 
temperature distribution, the observations of 
G333.6— 0.2 put a number of requirements on the 
parameters of any model. There must be a core 
of very high density to produce the bright radio 
and mid-IR core and the high [O III] 52/88 ^m 
and [S III] 19/33 /im line ratios seen in the center 
position. The material at 2' — 4' must be of low 
density to produce the 52/88 ^im and 19/33 /jm 
line ratios of order unity (N e of order 150 — 400 
cm -3 ). The material further out must be of even 
lower density (N e ~ 40 cm~ 3 ) to make a large 
contribution to the 122/205 /im line ratio ~ 2.5. 
The stellar T e g has to be high enough to produce 
++ and N ++ at all positions (out to > 3') but 
low enough so that there is ample N+, Ar+, and 
not much N ++ , ++ , or S 3+ even in the center 
position. In addition, the He + /H + radio recom- 
bination line (RRL) ratio is only 0.045 (McGee & 
Newton 1981), indicative of a relatively low T c tf 
(< 35 kK). 

We computed over a hundred models of G333.6— 0.2 
with a variety of combinations of stellar atmo- 
sphere models and density filling factors, all mod- 
els reproducing the observed radio flux. The abun- 
dance ratios used in the models were estimated 
from the fluxes in Table 4 (see also Table 3) and 
range from 1.25 to 3.4 times the Orion Nebula 
abundances (Simpson ct al. 1998). Line fluxes 
were computed from the models by integrating 
over flat-topped circular beam profiles for the IR 
lines or circular Gaussian beam profiles for the 
radio recombination lines located at the projected 
distances of the observed positions from the H II 
region center. 

Our best example of a spherically-symmetric 
model is the G333.6— 0.2 Core model in Table 4. 
This particular model is shown because it predicts 
line ratios that agree reasonably well with the den- 
sity indicators QS III] 19/33 ^m, [O III] 52/88 ^m, 
and [N II] 122/205 /xm) and the ionization indica- 
tor measured in the largest beam (He + /H + RRLs 
measured by McGee & Newton 1981). The model 
uses a linear combination of 0.9375 times the D-30 
(30 kK Dwarf) plus 0.0625 times the D-35 stellar 
atmosphere models of Pauldrach et al. (2001) and 
has a volume filling factor of 0.5. The contribution 
from the D-35 model, although seemingly small, 



is necessary to have enough high energy photons 
to ionize He and N to the observed ionization. 
The given model combination only represents the 
shape of the ionizing spectrum and not the actual 
numbers or spectral types of the ionizing stars. 

However, no model that we computed is a good 
fit — the data seem to require stellar ionizing 
fluxes with both higher and lower average T e g. 
The large fluxes of low-excitation ions, such as 
Ne + and Ar + in the center of G333.6— 0.2 seem 
to require stars with low T e ff (< 32 kK was sug- 
gested by Rank et al. 1978). However, the central 
density and optical depth are so high that when 
the exciting star is this cool, all the high energy 
photons that are capable of ionizing + to ++ 
are absorbed at small radii (0.87 pc, corresponding 
to 60") and the models have essentially no ++ 
at the radii corresponding to our three North posi- 
tions. Yet we observe ++ at all four positions. It 
would seem that an exciting star with a somewhat 
higher T c tf (> 35 kK) is required or a combination 
of stars including a much hotter star; however, 
for such models the > 35 eV photons required to 
ionize + to ++ still are confined in the cen- 
tral core but the total [O III] line fluxes are much 
larger than observed unless the O/H, O/N, O/Ne, 
etc. ratios are unreasonably low. On the other 
hand, in H II regions ionized by stars with rel- 
atively low T ff, the N ++ Stromgrcn radius, Rs, 
is much larger than the ++ Rs, leading to ob- 
served N ++ /0 ++ ratios that can be much larger 
than the actual N/O ratio (e.g., Rubin et al. 1994; 
SCREH; Stasihska & Schaerer 1997). The result 
of the differences in the N ++ and ++ Rs is that 
the observed N ++ /0 ++ ratio should increase with 
projected distance from the center. This predicted 
increase is not observed in our N ++ and ++ data 
taken at the North positions. Again, the observed 
behavior of this line ratio can be reproduced only 
by models with higher T e s stars. In fact, all the 
models that we computed predict a much larger 
[N III] 57/ [O III] 52 /im ratio than is observed at 
any position. Agreement cannot be improved by 
increasing the O/H ratio (the models already pre- 
dict more [O III] 52 /im flux than is observed) or 
decreasing the N/H ratio (the model [N II] lines 
fluxes agree with the observations). A lower T c s 
seems to be required, but already there is too lit- 
tle ++ in the North positions. We conclude that 
a spherical model, even of varying density, cannot 
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satisfy the observations of both the very bright 
core and the doubly ionized lines in the low den- 
sity region at 2' — 4' from the center. 

A possible solution might be some sort of non- 
spherical model with large clumps, or a blister 
model (face-on from its symmetrical appearance, 
Hyland et al. 1980). The NEBULA code for a 
blister model (Rubin et al. 1991) assumes that 
the density distribution is that of the surface of 
a molecular cloud with some possible gradient of 
density going into the cloud. The stellar photons 
ionize a hole on the edge of the cloud. However, 
that axisymmetric model does not use a density 
function that allows for a particularly high den- 
sity immediately next to the stars, as is required 
for a model of G333.6— 0.2. Without a modifica- 
tion to the axisymmetric 2-D code, this code is not 
suitable. A treatment that can handle a two- or 
three-dimensional density distribution appears to 
be required. 

We next considered a quasi-blistcr model or 
two-component model, such as was used by Simp- 
son et al. (1986), Colgan et al. (1991), and Moris- 
set et al. (2002) to model the Orion Nebula, K3- 
50, and G29.96-0.02, respectively. This would 
represent either a blister H II region with one open 
side or an H II region with very large clumps where 
ionizing photons escape through holes or paths be- 
tween the clumps. A linear combination was made 
of two spherical models, one with the high density 
core and one with an almost empty center to pro- 
vide the doubly ionized ions seen at 2' — 4'. The 
model with the empty center had a very large Rs 
and therefore had to be arbitrarily truncated in or- 
der that it be density bounded at about the same 
radius as the Rs for the Core model. An example 
of such a shell model, labeled "Model for the North 
Positions" , is also given in Table 4; the model used 
the same stellar atmosphere as the Core model but 
the filling factor was unity. Even here there are 
not enough photons to ionize ++ at the N2 and 
N3 positions and too much N ++ at the Nl posi- 
tion. Possibly three or even four spherical models 
with different interior shell diameters and differ- 
ent scaling fractions for the volume of the model 
in the line of sight are needed to reproduce the ob- 
servations at the four observed positions. Because 
of the arbitrariness of the procedure and the non- 
uniqueness of the result, we did not calculate any 
of these multi-shell models. 



The abundances have been only partially ad- 
justed to produce the observed fluxes — the imme- 
diate goal of the models was to determine T c g, the 
filling factor, and other geometry factors. Further 
adjustment to the abundances is not meaningful 
because the geometry still is not well determined, 
although the agreement of the He + /H + RRL ra- 
tios with the observations indicates that T e g is not 
especially discrepant (since most helium is primor- 
dial, the interstellar He/H ratio has not increased 
significantly as a result of Galactic chemical evo- 
lution). The uncertainty in the geometry factors 
makes it difficult to estimate a reliable abundance 
for oxygen. However, if we consider just the Cen- 
ter and Nl positions, we find that the observed 
fluxes would be matched by the model fluxes if 
O/H^ 3 x 10~ 4 , a value somewhat lower than that 
observed in the Orion Nebula. This O/H abun- 
dance ratio is not reliable, but suggests that O/H 
is not especially high in G333.6— 0.2. As expected, 
since we observed the dominant ionization states 
of nitrogen and sulfur (N + and S ++ ) in this source, 
our final estimated abundances of these elements 
are those input to the models (Table 3). On the 
other hand, it is difficult to think of a geometry 
that would produce the large [Ar II] 6.98 /im line 
flux observed in the G333.6— 0.2 core, given the 
number of high energy photons that arc needed to 
produce the observed [Ne III] and [O III] lines. 
This discrepancy may indicate that the atomic 
data for argon need revision (see also Morisset et 
al. 2004). 

4.2.2. W43 

Because the constant density model by SCREH 
for W43 is not a good match to the densities in- 
ferred from the observations (see Table 2), we cal- 
culated new models for it. The appearance of W43 
is that of a large "T" or incomplete "D" shape 
projected on the sky (Lester et al. 1985; Liszt 
1995; Two Micron All Sky Survey image 5 ) with 
a central star cluster near but not on the vertical 
element of the letter. Because our three observed 
positions, including the one closest to the cluster, 
all have similar densities and excitation (SCREH; 
Table 2), we interpret the geometry as a cluster of 
ionizing stars surrounded by a cavity and then a 
shell of ionized gas, such that all three positions 

5 2MASS, http://irsa.ipac.caltech.edu 
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are at essentially the same physical distance from 
the stars. Thus we model the H II region as a shell 
of ionized gas with interior radius = 1.2 pc (0.75' 
at a distance of 5.6 kpc (Reifenstein et al. 1970; 
Liszt 1995), scaled to a Galactic Center distance 
of 8 kpc). The model N e drops smoothly from 
770 cm" 3 at 1.2 pc to 420 cm" 3 at R s ~ 2.9 pc, 
matching the average of the total observed flux of 
77.5 Jy at 2 cm (Schraml & Mezger 1969) and 97.4 
Jy at 6 cm (Reifenstein ct al. 1970) for an ionizing 
luminosity of 2.9 x 10 50 photons s _1 . The densities 
and ionic abundance ratios of the three positions 
were averaged (weighted) for comparison with the 
models; these comparisons are given in Table 5. 

Three models were computed with the only dif- 
ference being the input stellar atmosphere mod- 
els. The stellar atmosphere used for Model D is 
the sum of 0.8 times the D-35 model atmosphere 
from Pauldrach et al. (2001) plus 0.2 times the D- 
40 model atmosphere, the stellar atmosphere used 
for Model S is the sum of 0.6 times the S-35 model 
atmosphere plus 0.4 times the S-40 model atmo- 
sphere, where D-35, D-40, S-35, and S-40 stand 
for Dwarf and Supcrgiant models with T e g equal 
to 35 kK and 40 kK, respectively, and the stel- 
lar atmosphere used for Model T is the TLUSTY 
model atmosphere from Lanz & Hubeny (2003) 
for T e ff = 35 kK and log g = 3.50. The reason the 
S (supergiant) model atmospheres were also used 
is the presence in W43 of a Wolf-Rayet star and 
two O giants/supergiants (Blum et al. 1999). The 
combinations of stellar atmosphere models were 
adjusted until the H II region models produced 
the observed average radio He + /H + ratio (He/H 
= 0.10 for the nebular models) and the observed 
N++ /N+ ratio. 

To correct the observed ionic abundance ratios 
for unseen ionization states, we used the ioniza- 
tion fractions calculated from the model as shown 
in the Table 5 comments and footnotes: the aver- 
aged ionic abundance ratios were divided by the 
model ionic ratios, which are inverse ICFs. For 
example, the average observed ++ /S ++ ratio is 
16.33 and Model S has <0++/0>/<S++/S>= 
0.391. Using this inverse ICF, we get O/S = 
41.8 and O/H = O/S * S/H = 4.8 x 10~ 4 , using 
S/H= 11.4 x 10~ 6 for Model S from the average 
S++/H+ = 10.9 x 10~ 6 . The main result to be 
noticed is that the N/H ratios are fairly similar 
for all models at 1.2 to 1.5 x 10~ 4 , even though 



the O/H ratio varies by a factor of 2.5. The N/S 
ratio calculated by averaging N/S from iV + /S ++ 
and N++/S++ is 9.9 to 10.4, very similar to the 
result from summing (N ++ + N + )/S ++ found in 
Table 3. 

The O/H variation is caused by the model 
<0 ++ /0> ionization fractions, as seen in Table 5 
in the <0++/0>/<S++/S> and <N++/N>/<0++/0> 
ratios. No doubt the reason is that the Supergiant 
model atmospheres have less flux from 35 to 41 eV 
than the Dwarf model atmospheres. Pauldrach et 
al.'s (2001) Supergiant models have much higher 
mass loss rates than their Dwarf models, so the 
difference must be caused by more line blanketing 
at these energies by the gas in the wind. On the 
other hand, the static, plane-parallel, non-LTE at- 
mospheres of Lanz & Hubeny (2003), which have 
no winds, have much lower fluxes from ~ 40 to 
- 50 eV. 

All three W43 models adequately reproduce the 
T cff -indicating indices: the He+/H+ and N++/N+ 
ratios. The range of implied abundances in Ta- 
ble 5, though, indicates the large uncertainties re- 
sulting from determination of abundances through 
model-produced ICFs, particularly when different 
stellar atmosphere models are used to compute the 
ICFs. We prefer the values of N/H derived from 
comparisons with sulfur, rather than those derived 
from oxygen. The estimated N/O ratio is 0.32 
to 0.42 from ICFs computed with models D and 
S (not very different from the 0.34 estimated by 
SCREH), but is as low as 0.16, similar to that of 
the Orion Nebula (Rubin et al. 1991), from the 
ICFs computed with model T. The reason for the 
wide range is not the abundances of nitrogen or 
sulfur, which are reliable because they do not re- 
quire large ICFs. The problem is the total abun- 
dance of oxygen: the model ICFs produce O/H 
ranging from 2.8 to 8.3 xl0~ 4 , whereas the pre- 
dicted O/H ratio is 7.2 x 10 -4 for a Galactic abun- 
dance gradient of d(\og 0/H)/dR — —0.06 dex 
kpc -1 (Henry & Worthy 1999), given the Orion 
Nebula's location at R G ~ 8.4 kpc and O/H = 4.0 
to 4.4 x 10~ 4 (Rubin et al. 1991; Esteban et al. 
1998). In fact, one expects that the oxygen abun- 
dance for W43 should be much higher than that 
of the Orion Nebula, given W43's low T e = 5410 
to 6500 K (Subrahmanyan & Goss 1996; Wink, 
Wilson, & Bieging 1983; and references in both 
papers). For this source, as for G333.6— 0.2, we 
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conclude that we still do not have a good estimate 
for the abundance of oxygen, the most abundant 
heavy element and the most important contribu- 
tor to cooling in H II regions. 

5. Discussion and Conclusions 
5.1. Stellar Atmosphere Models 

As was discussed above, many of the ioniza- 
tion correction factors are critically dependent on 
the choice of stellar atmosphere model. This can 
be demonstrated either by plotting ratios of the 
numbers of photons able to ionize various ions as 
a function of T e g- or by plotting ionization frac- 
tions or ICFs calculated from models; here we do 
the latter because ICFs are more readily appli- 
cable to abundance computation. This is illus- 
trated in Figure 5, where we plot the ICF ratios 
<N++/N>/<0++/0>, <0++/0>/<S++/S>, 
and <0++/0> versus the <N++/N>/<N+/N> 
ratio, which indicates excitation without the ne- 
cessity for any abundance assumption. The mod- 
els from Rubin (1985) that are plotted in Fig. 5 
all have total nucleon density = 10 3 cm~ 3 , ion- 
izing luminosities of 10 49 to 10 50 photons s _1 , 
and Orion Nebula metallicities; the models from 
SCREH had similar densities and ionizing lumi- 
nosities but metallicities ranging from 1 to 2 times 
Orion nebula metallicities. The models from both 
Rubin (1985) and SCREH were computed with 
the versions of NEBULA available at the time 
and both used LTE stellar atmosphere models 
from Kurucz (1979). In addition to the models 
described in §4, constant nucleon density (= 10 3 
cm~ 3 ) spherical models with ionizing luminosities 
of 10 49 photons s _1 were computed and plotted. 
These recent NEBULA models used stellar atmo- 
sphere models from either Pauldrach et al. (2001), 
Lanz & Hubcny (2003), or contemporary 6 LTE 
stellar atmosphere models described by Kurucz 
(1992). The stellar atmosphere models had "So- 
lar" abundances and the nebular models had 1.5 
times Orion Nebula abundances (Simpson et al. 
1998) such that the model abundance ratios (all 
times 10~ 4 ) are C/H = 3.75, N/H = 1.02, O/H 
= 6.0, Nc/H = 1.5, S/H = 0.11, and Ar/H = 
0.04. Because they are widely used, we also plot 
ICFs from the models of Stasinska & Schaerer 



6 http: / /kurucz. harvard.edu/ 



(1997). For these models both the stellar and 
nebular abundances were "Solar". We note that 
some of the elemental abundances considered to be 
"Solar" are changing with time; for example the 
"Solar" O/H ratio used by Stasinska & Schaerer 
(1997) was 0/H= 8.51 x 10~ 4 (Anders & Grevesse 
1989), whereas currently, the "Solar" O/H ratio 
is thought to be 5.45 x 10~ 4 (Holweger 2001). 
Martin-Hernandez et al. (2002b), Smith et al. 
(2002), and Morisset (2004) discuss the effects of 
stellar abundances on H II ionization. Details of 
the physics used, though, are much more impor- 
tant than the exact "Solar" abundances, as we 
will see. 

Unless an H II region can be observed at both 
optical and infrared wavelengths, observers cannot 
measure the abundances of N, O, Nc, S, and Ar 
without making corrections for unseen ionization 
states for some elements. For optically obscured 
H II regions like W43 and G333.6-0.2, the most 
significant missing ion is + , since oxygen is the 
most abundant heavy element and since these low 
excitation H II regions are believed to have most 
of their oxygen as + . Rubin et al. (1994) and 
SCREH assumed that they knew the O/S ratio 
and used the observed ++ /S ++ ratio to estimate 
the ICF for their measured N ++ /0 ++ ratios, since 
they did not have any other line pair that could 
be used to indicate excitation. In this paper we 
use our observed N ++ /N + ratios to estimate T e g; 
however, we see from Fig. 5a how dependent the 
N ++ /0 ++ ICF is on the stellar atmosphere model. 

Martin-Hernandez et al. (2002a) took ad- 
vantage of their ISO SWS measurements of the 
Ne ++ /Ne + ratios to estimate ICFs without the 
need for assuming a priori some abundance ratio; 
unfortunately, the large beam size difference be- 
tween the SWS and LWS instruments means that 
their observed Ne ++ /Ne + ratios may not always 
be appropriate for estimating the excitation of 
++ and N ++ . Moreover, Ne ++ requires ionizing 
photons > 41 eV, and it is at these high energies 
that the stellar atmosphere models are especially 
divergent, as we now discuss. 

An important test of the validity of stellar at- 
mosphere models of hot stars is whether H II re- 
gion models produced with these atmospheres pre- 
dict line fluxes that agree with observations. Ex- 
cept for a few rare instances, this is the only way 
that the EUV fluxes of the models can be tested, 



12 



since interstellar absorption by hydrogen predom- 
inantly prevents the EUV fluxes from hot O stars 
from reaching the Earth. Here we discuss the stel- 
lar fluxes that can doubly ionize neon (photon en- 
ergy > 41 cV). This has been called the [Ne III] 
problem (e.g., Sellmaier et al. 1996) because it has 
been observed that the Ne ++ /0 ++ ratio is rela- 
tively constant over a large range of H II region 
excitation (SCREH; Stasihska & Schaerer 1997 
and references therein; Kennicutt et al. 2003), 
contrary to predictions of models that used LTE 
stellar atmospheres (e.g., SCREH). Proposed so- 
lutions have involved using non-LTE stellar atmo- 
spheres, usually with winds (Rubin, Kunze, & Ya- 
mamoto 1995; Sellmaier et al. 1996; Stasihska 
& Schaerer 1997). Figure 6 shows the observa- 
tions of SCREH along with the predictions of the 
models used in Fig. 5 and the G333.6— 0.2 and 
W43 models described in §4. Models calculated 
with black bodies for the stellar atmosphere spec- 
trum are also plotted since Morisset et al. (2004) 
found that such models give a reasonable fit to 
the mid-IR ISO observations. Like the new neb- 
ular models plotted in Fig. 5, these models also 
have spherical geometry, constant nucleon density 
of 10 3 cm -3 , 10 49 ionizing photons, and 1.5 times 
the Orion Nebula abundances listed in Table 3. 

Surprisingly, in Fig. 6 only the models of 
Stasihska & Schaerer (1997) reproduce the obser- 
vations over the total range of observed ++ /S ++ , 
although if Nc/O is as large as 0.25 in all H II re- 
gions as it is in the Orion Nebula (Simpson et al. 
1998; Table 3), these H II region models produce 
too high a Ne ++ /0 ++ ratio because of the very 
large EUV fluxes in Schaerer & dc Kotcr's (1997) 
stellar atmosphere models. Moreover, Schaerer & 
de Koter's models have too much EUV flux for 
reproducing other observations (e.g., Smith et al. 
2002). Earlier, Sellmaier et al. (1996) had demon- 
strated that the [Ne III] problem was supposedly 
solved when they obtained a good fit to the data 
by using non-LTE atmospheres with winds com- 
puted with Pauldrach's code as it then existed. 
However, our H II region models with non-LTE 
stellar atmospheres with winds from Pauldrach et 
al. (2001) taken from both A. Pauldrach's web site 
and from the stellar atmosphere models produced 
by Smith et al. (2002) using Pauldrach's WM- 
basic code predict much lower Ne ++ /0 ++ than 
observed for Dwarf atmospheres with T c s < 40 



kK and for Supergiant atmospheres with T ff < 35 
kK. Our H II region model produced with a 35 kK 
supergiant atmosphere from F. Sellmaier (private 
communication in 1995) lies above the Pauldrach 
et al. Supergiant line but not as high as the model 
with the 35 kK supergiant atmosphere in Sellmaier 
et al.'s (1996) Figure 2, even though we are now as- 
suming Ne/O = 0.25 (Simpson et al. 1998) instead 
of 0.2025. Part of the difference of the current neb- 
ular model predictions with previous models is the 
higher ++ ionization relative to Ne ++ and S ++ 
obtained with the Opacity Project cross sections 
described in Section 4. The model computed with 
Scllmaicr's (1995) atmosphere has higher ioniza- 
tion than the model computed with Pauldrach et 
al.'s (2001) S-35 atmosphere because it used Orion 
Nebula abundances (as in Sellmaier et al. 1996) in- 
stead of the higher mctallicity (1.5 times the Orion 
Nebula abundances of Table 3) used for the rest of 
the models plotted here. The result, though, is a 
return of the mismatch between the stellar atmo- 
sphere models and the EUV fluxes of real stars, 
with possible additional discrepancies for the ion- 
ization, recombination, and collisional excitation 
cross sections for oxygen, neon, and sulfur. More- 
over, we see that all the models computed with- 
out winds do a poorer job of reproducing the ob- 
servations than both the models computed with 
non-LTE atmospheres with winds and the models 
computed with simple black bodies for stellar at- 
mospheres (for this reason we give low weight to 
the T model for W43 in Table 5). This probably 
means that the codes producing the stellar atmo- 
sphere models used for computation of H II region 
models should include the effects of stellar winds. 

5.2. Abundance Ratios 

At the low excitations seen in inner Galaxy 
H II regions, the predicted flux of the ++ lines 
is extremely sensitive to the stellar atmosphere, 
but if the models are accurate, the O/H ratio of 
both G333.6— 0.2 and W43 could be as low as 
3 — 5 x 10 -4 , substantially lower than would be 
predicted from the observed Orion Nebula abun- 
dance and an O/H abundance gradient of —0.06 
dex kpc" 1 (Henry & Worthey 1999). This low 
O /H ratio is hard to understand because of the low 
T e measured in inner Galaxy H II regions from ra- 
dio recombination lines. There would need to be 
some source of systematic error for the O/H ra- 
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tio to be higher than these estimates: i.e., what 
would be required is that the excitation of oxy- 
gen in the two H II regions is lower than esti- 
mated in this paper from He + /H + recombination 
lines and N + and N ++ forbidden lines. Possible 
reasons that the G333.6— 0.2 and W43 excitation 
might be overestimated are as follows: (1) The 
abundance of N + is underestimated because the 
[N II] line fluxes are undermeasured owing to tel- 
luric absorption at 122 /im or diffraction at 205 
fim. This is a particular problem for W43, where 
the two measurements of the [N II] 122 /im line in 
different apertures can be used to produce very dif- 
ferent electron densities and N++/N+ abundance 
ratios (0.6 to 3.2 in Table 5). (2) The estimated 
N ++ /N + ratio is too high because some of the 
estimated N + is missing owing to inadequate cor- 
rection for density variations. (3) The estimated 
N ++ /N + ratio is too high because of incorrect col- 
lisional excitation cross sections for the pertinent 
energy levels. (4) The models predict too much 
0++/0 versus N++/N+. The spread of 0++/0 
ratios in Fig. 5c shows that any excessive ++ /0 
is due to model atmospheres and not due to the 
~ 30% higher oxygen ionization resulting from the 
Opacity Project cross sections. (5) The ioniza- 
tion of helium might be overestimated if the Hc/H 
abundance ratio is actually greater than the as- 
sumed 0.10. We conclude that determinations of 
the O/H abundance ratio in low excitation H II 
regions from FIR data have large uncertainties at 
this time, with the largest contributor to the un- 
certainty being the choice of stellar atmosphere 
models, with some possible contribution from the 
atomic data used in the nebular models. 

We have demonstrated that the abundances of 
nitrogen and sulfur can be obtained from FIR 
observations and the ratio of N/S can be ob- 
tained when both ionization states of N + and N ++ 
are measured. (Simpson et al. 1998, Martin- 
Hernandez et al. 2002a and 2003, and Giveon et 
al. 2002 also show that the abundances of Ne/H 
and Ar/H, as well as S/H, can be determined from 
mid-IR observations with little bias owing to exci- 
tation or abundance.) Compared to SCREH, the 
S/H abundance ratio is higher in G333.6— 0.2, ow- 
ing to the new data at the additional North po- 
sitions, but lower in W43 because of the higher 
collisional excitation cross sections for sulfur. The 
N/H ratio is about 50% larger in G333.6— 0.2, ow- 



ing to the new data at the North positions, but a 
factor of 2 smaller in W43, owing to the smaller 
ICF for N + . However, since the estimated O/H ra- 
tio is also much smaller in W43, now that we are 
estimating the ionization from N ++ /N + and not 
from ++ /S ++ with the assumption of a constant 
O/S ratio as SCREH did, the estimated N/O ra- 
tios agree with SCREH for both G333. 6-0.2 and 
W43. With this new analysis, the S/H and N/O 
ratios are also consistent with SCREH's measure- 
ments. 

5.3. Conclusions Regarding Systematic 
Errors 

There are several possible sources of system- 
atic error that must be considered when measur- 
ing abundances in obscured, low-excitation H II 
regions, that is, all the H II regions in the inner 
Galaxy. The impact of these systematic errors is 
that there is a very large uncertainty in the total 
metallicity and important abundance ratios like 
N/O. These sources of systematic error are the 
following: 

(1) The most abundant heavy element ion in 
low-excitation H II regions, + , has no bright in- 
frared lines. The consequence is that the abun- 
dance of this critically important clement must 
be estimated from the abundance of an ionization 
state with fractional abundance < 0.5. 

(2) Nitrogen, the secondary element which one 
needs for studies of galaxy chemical evolution, has 
two important ionization states: N + as well as 
N ++ . Fortunately, both species have measurable 
FIR lines, but both must be measured to determine 
the excitation. 

(3) Moreover, both FIR N + lines have greatly 
different values of N cr n from those of S ++ or N ++ . 
Thus one needs to be aware of and compensate for 
the ever-present density variations. 

(4) The average density could be quite differ- 
ent in the various ionization zones, as shown by 
the different electron densities derived from N e - 
diagnostic line pairs: low ionization (N + ), in- 
termediate ionization (S ++ ), and high ionization 
(0++). 

(5) ICFs estimated from models have major 
uncertainties owing to uncertainties in the stellar 
atmosphere models and possibly the atomic data 
used in the H II region models. 
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It is clear that the abundances of primary Astrophysics grant, which funds his Co-operative 

(oxygen, neon, or sulfur) and secondary (nitro- Agreement NCC2-9021. 

gen) elements will need to be measured in many 
more inner Galaxy H II regions before we under- 
stand the chemical evolution of the Milky Way. 
In the near future, measurements of N ++ , N + , 
++ , Nc ++ , and S ++ will be possible from the 
Stratospheric Observatory for Infrared Astron- 
omy, SOFIA. We recommend that future observers 
measure excitation-sensitive lines, such as Ne + 
and Ne ++ , as well as both N+ and N ++ to de- 
termine the H II region ionization state needed to 
estimate the ICF for + . It is important that a 
sufficient number of lines be measured in each H II 
region (and extended H II regions mapped) so that 
detailed models of the H II region can be made. 
Discrepancies between the predicted line fluxes 
from the models and the observations can then 
be used to indicate EUV energy regions where the 
stellar atmosphere models may need revision. 
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Fig. 1 . — (a) Ratios of volume emissivities for line 
pairs of the electron-density-sensitive lines [S III] 
19/33 fiia, [O III] 52/88 fj,m, and [N II] 122/205 
/zm plotted as a function of electron density, N e , 
for electron temperature, T e , of 7000 K. (b) Line 
intensities were calculated for slab models with 
7V e = Noe~ x by integrating the volume emissiv- 
ity of each line from Xq to x max (see text). Values 
of N e were estimated from the abscissa of Fig. la 
by assuming that the ratios of the line intensities 
could be used for the ordinate. These estimated 
values of N e divided by No are plotted versus No. 
The root-mean-squared (RMS) density divided by 
TVo is also plotted (dotted line). 
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Fig. 2. — Normalized abundance measures calcu- 
lated for slab models with exponentially decreas- 
ing density (see text) plotted as a function of the 
densities estimated from the ratios of the density- 
sensitive lines (see Fig. la). The abundance mea- 
sures, AM em and AMcd, are defined as the cal- 
culated line intensity divided by the normalized 
volume emissivity or the emissivity per ion, re- 
spectively, and are normalized by dividing by the 
emission measure or by the column density of ions, 
respectively (see text). Because preferred abun- 
dance measures for each ion depend on the den- 
sity, critical densities, N cr u, are also tabulated for 
T e = 7000 K. (a) S++ measured from [S III] 19 
and 33 /im lines, (b) ++ and N ++ measured from 
[O III] 52, [O III] 88, and [N III] 57 lines. Note 
the near coincidence of the N ++ and ++ abun- 
dance measures, (c) N + measured from [N II] 122 
and 205 /tm lines. 




Fig. 3. — Ratios of the abundance measures for the 
ions shown in Fig. 2 plotted as a function of esti- 
mated N e for lines of sight with varying amounts 
of high-density gas. The heavy line is AMem and 
the dotted line is AMcd (the light line is unity, for 
reference), (a) N++/0++. The lines for AM EM 
and AM CD coincide, (b) 0++/S++. (c) N+/S++. 
(d) N++/S++. (e) N++/N+. 



21 




Fig. 4.— (a) Azimuthally-averaged radio bright- 
ness temperatures for G333.6— 0.2 plotted as a 
function of angular distance on the sky from the 
central peak. The radio brightness temperatures 
were all scaled to the same frequency, 5 GHz, us- 
ing the formulation of Mezger & Henderson (1967) 
and ignoring optical depth effects, (b) N e , derived 
by inverting the radio brightness temperature dis- 
tribution shown in Fig. 4a, is plotted as a func- 
tion of distance, r, from the central star cluster. 
A spherical H II region at a distance of 3 kpc was 
assumed. The overall form of this density function 
is N e oc r -12 , as shown by the dashed line, which 
was fit by least squares. 



22 



Fig. 5. — Ratios of ionization fractions plotted 
for spherical H II region models. The models 
with LTE stellar atmosphere models are those 
from Rubin (1985) and SCREH, who used model 
atmospheres from Kurucz (1979), and the new 
NEBULA models with contemporary atmosphere 
models described by Kurucz (1992) (T c s = 33, 
35, 37, 40, and 45 kK). The H II region models 
with non-LTE stellar atmosphere models are those 
from Stasihska & Schaerer (1997), who used atmo- 
sphere models from Schaerer & de Koter (1997), 
and the NEBULA models calculated for this pa- 
per, where we used non-LTE atmosphere mod- 
els from Pauldrach et al. (2001) and Lanz & 
Hubeny (2003) (T eff = 30, 35, 40, and 45 kK). 
For the H II region models calculated with Paul- 
drach ct al. atmospheres, the solid line con- 
nects models with Dwarf atmospheres and the 
dashed line connects models with Supcrgiant at- 
mospheres. For the H II region models cal- 
culated with Lanz & Hubeny atmospheres, the 
solid line connects models with atmospheres with 
log g = 4.0 and the dotted and dashed lines 
connect models with atmospheres with log g = 
3.0 to 3.5 and with Lyman continuum luminosi- 
ties of 10 49 and 10 50 photons s _1 , respectively. 
The notation in the legend applies to all parts 
of the figure, (a) <N++/N>/<0 ++ /0> versus 
<N++/N>/<N+/N>. (6) <0++/0>/<S++/S> 
versus <N++/N>/<N+/N>. (c) <0++/0> ver- 
sus <N++/N>/<N+/N>. 




0.01 0.10 1.00 10.00 100.00 

<N*7N>/<N7N> 
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Fig. 6. — Ratios of Ne ++ /0 ++ inferred from 
the observations of Simpson et al. (1995) plot- 
ted against their inferred ++ /S ++ ratios. Ra- 
tios from the nebular models discussed in this pa- 
per are also plotted. The notations of Pauldrach 
ct al. (2001), Sellmaier (private communication, 
1995), Lanz & Hubeny (2003), and Kurucz (1992) 
in the legend refer to the stellar atmospheres used 
in our NEBULA models (same T c ff as the mod- 
els in Fig. 5; Sellmaier's atmosphere is a Super- 
giant with T cS = 35 kK). The plotted Ne++/0 ++ 
and ++ /S ++ ratios for the models are equal 
to the calculated <Ne ++ /Ne>/<0++/0> and 
<0++/0>/<S ++ /S> ratios multiplied by as- 
sumed abundance ratios of Ne/O = 0.25 and O/S 
= 50, respectively (Orion Nebula abundance ra- 
tios from Simpson et al. 1998). For the H II mod- 
els calculated with Pauldrach et al. atmospheres, 
the solid line connects models with Dwarf atmo- 
spheres and the dashed line connects models with 
Supcrgiant atmospheres. For the H II region mod- 
els calculated with Lanz & Hubeny atmospheres, 
the solid line connects models with atmospheres 
with log g = 4.0 and the dotted and dashed lines 
connect models with atmospheres with log g = 3.0 
to 3.5 and with Lyman continuum luminosities of 
10 49 and 10 50 photons s _1 , respectively. Models 
computed with atmospheres represented by black 
bodies from 20 kK to 40 kK (every 5 kK) are also 
plotted (open circles connected by a dotted line). 
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Table 1 

Observed Line and Continuum Fluxes. 



Lino Flux" (10~ 18 W cm" 2 ) 
Continuum Flux" (Jy) 
G333.6-0.2 W43 (G30. 8-0.0) 



Line 


G333.6-ISO b 


G333.6-N1 C 


G333.6-N2 C 


G333.6-N3° 


W43N C 


W43C C 


W43S C 


[S III] 18.7 




160.1±16.2 


36.3±6.4 




SCREH 


SCREH 


SCREH 




1370±275 


290±280 




SCREH 


SCREH 


SCREH 


[S III] 33.5 (im 




189.1±3.6 


65.7±4.7 


38.9±3.1 


SCREH 


SCREH 


SCREH 






7040±580 


740±185 


730±180 


SCREH 


SCREH 


SCREH 


[O III] 51.8 (im 


178.2±8.5 


84.3±2.6 


13.1±1.6 


6.8±1.2 


SCREH 


SCREH 


SCREH 




54400 


9910±320 


1130±360 


410±270 


SCREH 


SCREH 


SCREH 


[N III] 57.3 Aim 


46.3±4.9 


51.1±2.8 


11.9±1.4 




SCREH 


SCREH 


SCREH 


57630 


10945±230 


1310±200 




SCREH 


SCREH 


SCREH 


[O III] 88.4 


39.9±3.9 


37.7±0.8 


14.5±0.9 


7.5±0.7 


SCREH 


SCREH 


SCREH 


55100 


9030±175 


1060±270 


630±90 


SCREH 


SCREH 


SCREH 


[N II] 121.9 nm 




7.4±1.4 


2.6±0.8 




14.4±0.8 


5.5±0.5 


7.9±0.8 




6790±120 


575±230 




5270±420 


1450±260 


4690±380 


[N II] 121.9 ^im d 










4.9±0.5 d 
2530±210 d 


2.5±0.4 d 
800±200 d 


2.6±0.6 d 
1600±300 d 


[N II] 205.2 ^irn 




2.0±0.2 


1.5±0.2 


1.7±0.2 


1.52±0.26 


0.46±0.14 


0.56±0.13 




2000±145 


230±120 


115±100 


1300±70 


490±50 


1440±40 



a Flux and continuum errors arc 1 a and are statistical only. 
b All ISO LWS data had entrance aperture 80". 
C A11 CGS data had entrance aperture 60" except as marked. 
d Entrance aperture 45". 
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Table 2 

Electron Densities and Abundance Ratios 







G333.6-0.2 








W43 




G333.6-ISO 


G333.6-C 


G333.6-N1 


G333.6-N2 


G333.6-N3 


W43N 


W43C 


W43S 


T9.6 1.5 
S 5G Hz b 40.75 


1.5 
26.65 


1.5 

5.35 


1.5 
1.10 


1.5 
0.35 


2.74 
5.91 


2.74 
2.25 


3.63 
3.11 


Electron Densities (cm - 3 ) 
JV e (S III) 

N e (O III) 2520±540 e 
N e (N II) 

N (N TT^l . . . 


2400±350 c 
4200±700 c 

,, n + 475c 


735±155 d 
730±50 d 

, 4r) + 6Sd 


335±150 d 
140±50 d 

oy+24d 
31 -21 


140±70 d 


340±120° 
870±50 c 

5120^°° d 
1 enxenf 

IOUIOU 


380±130 c 
730±60° 

oo d 


280±160 c 
725±50 c 
oo d 

42n +835f 


Ionic Abundance Ratios (xlO -6 ) 
S + + /H+ 

0++/H+ 38±5 
N++/H+ 14±2 


4.3±0.7 
71±9 
36±4 


12.2±1.6 

80±8 
48±5 


14.2±2.1 
50±6 
27±4 


20.8±2.7 
82±12 


9.0±1.4 
130±14 
68±8 


14.1±2.2 
275±31 
127±14 


12.8±2.3 
227±25 
117±13 


Ionic Abundance Ratios 

N+ + /O++ 0.36±0.04 
O++/S++ 
N++/S++ 
N+/S++ 
N++/N+ 


0.51±0.02 
16.7±2.2 
8.47±0.70 
0.63±0.20 
13.3±4.0 


0.60±0.04 
6.53±0.55 
3.90±0.29 
1.90±0.52 
2.06±0.45 


0.53±0.07 
3.56±0.57 
1.89±0.29 
1.75±0.45 
1.08±0.27 


3.95±0.52 
14.3±1.3 


0.52±0.03 
14.5±1.2 
7.54±0.38 
1.83±0.37 g 
4.11±0.72 g 


0.46±0.02 
19.5±1.9 
9.04±0.36 
3.98±1.67 g 
2.27±0.84 g 


0.52±0.02 
17.8±1.8 
9.18±0.37 
2.30±1.38 g 
3.99±1.46 g 


Column Density Ratios 

O++/S++ 
N++/S++ 

N+/S++ 

N++/N+ 

N+/S++ 

N++/N+ 


9.55±0.43 
4.85±0.23 
4.62±0.58 
1.05±0.13 


6.56±0.22 
3.92±0.24 
9.96±1.29 
0.39±0.08 


8.60±0.87 
4.56±0.67 
15.7±2.4 
0.29±0.09 


3.95±0.52 
14.3±1.3 


5.64±0.22 

2.94±0.15 

4.61±0.36 h 

0.64±0.05 h 

3.39±0.41 g 

0.87±0.12 g 


10.2±0.5 
4.72±0.18 
3.29±0.37 h 
1.44±0.14 h 
2.31±0.37 g 
2.04±0.34 g 


6.89±0.47 

3.56±0.19 

1.95±0.50 h 

1.27±0.14 h 

1.55±0.34 g 

2.30±0.52 g 


Adopted Abundance Ratios 

O++/S++ 
N++/S++ 
N+/S++ 
N++/N+ 
(N++N+ + J/S++ 


13.1±3.6' 
6.66±1.8l' 
2.63±1.99' 
1.05±0.13 k 

9.3±3.8 


6.53±0.55 
3.90±0.29 
5.93±4.03' 
1.23±0.83' 
9.8±4.3 


3.56±0.57 
1.89±0.29 
8.71±6.96' 
0.69±0.39' 
10.6±7.3 


3.95±0.52 
14.3±6.0' 
>14±6 


14.5±1.2 
7.54±0.38 
2.96±0.38 j 
1.43±0.30' 

10.5±0.8 


19.5±1.9 
9.04±0.36 
3.02±0.80 j 
1.44±0.33' 

12.1±1.2 


17.8±1.8 
9.18±0.37 
2.59±0.74 j 

I. 89±0.53' 

II. 8±1.1 



a Ccnter position, data from Colgan ct al. (1993). 

b Radio flux in the aperture used for the FIR line measurement, scaled to the value appropriate for v = 5 GHz. 
c Elcctron density computed using the 45 /; aperture data. 
d Elcctron density computed using the 60" aperture data. 
e Electron density computed using the 80" aperture data. 

f Electron density computed using the 45 // aperture data for the 122 fim line and reducing the 205 fim line flux by the ratio of the 
aperture diameters. 

s Computed using the 45 aperture data for the 122 fim line and reducing the 205 fim line flux by the ratio of the aperture diameters. 

Computed using the 60" aperture data for the 122 /^m line and reducing both the 122 fim and the 205 fim line fluxes by the ratio of 
the aperture diameters for comparison to the [N III] and [S III] lines. 

1 Average of column density ratio and ionic abundance ratio; the error is increased for G333.6-N3 where N e (N II) was not measured. 

J Average of ionic abundance ratios for all 122 /j.m line measurements plus column density for the 45" aperture 122 /j-m line measure- 
ments. 

k Column density ratio chosen because of high N e compared to N crit . 
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Table 3 
Estimated Abundances 



Element 


G333.6-0.2 


W43 


Orion Nebula 


Estimated S/H a 


13.4±0.8 x 10~ 6 


i2.i±i.i x icr 6 


7.0 x 1CT 6 b 


Estimated N/S a 


9.3±2.1 


10.3±0.4 


9.71 b 


Estimated N/H c 


12.5xl(T 5 


12.5xl0~ 5 


6.8xl0~ 5 d 


Estimated Nc/H 


3.4xlO~ 4 


l.lxlO~ 4 f 


1.0xl0~ 4 b 


Estimated 0/H 


3 to 4xl(T 4 g 


3 to 8xlCT 4 s 


4.0xl(T 4 d 



a Assumed fractional ionization <S ++ /S>= 0.9. 

b Simpson et al. (1998) revised for new collisional excitation cross sec- 
tions (see text). 

C N/H = N/S x S/H. 

d Rubin et al. (1991). 

c Estimated from [Ne II] 12.8 /im measurements of Fujiyoshi et al. (1998) 
and radio continuum measurements of Fujiyoshi (1999). 

f SCREH. 

SO/H = 0++/H+ x ICF Modeh 
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Table 4 

G333.6-0.2 - Comparison of Models to Observations 



Line Flux or Ratio 


Observed Fluxes and Ratios" 


Model for G333.6-0.2 C 




Model for North Positions 11 ' 


Center 


Nl 


N2 


N3 


Center 


Nl 


N2 


N3 


Nl 


N2 


N3 


[O III] F 52 fim 


161° 


92 


14 


7 


350 


0.0005 








117 


1 





[O III] F 88 nm 


27 c 


39 


15 


8 


54 


0.0001 








67 


1 





[O III] F 5 2/F 8S 


6.06 


2.37 


0.96 


0.96 


6.50 


4.12 






1.75 


1.61 




[S III] F 19 fim 


412 c 


270 


61 




869 


273 


50 


14 


154 


52 


16 


[S III] F 33 


158 c 


228 


79 


47 


351 


243 


72 


27 


192 


85 


33 


[S III] F 19 /F 33 


2.61 


1.19 


0.78 




2.48 


1.13 


0.70 


0.52 


0.80 


0.61 


0.48 


[N II] Fi 22 Aim 


2.5 C 


7.4 


2.6 




9.2 


13.7 


5.6 


2.8 


6.8 


6.9 


3.5 


[N II] F 205 Mm 


0.44 c 


1.97 


1.46 


1.66 


1.5 


2.4 


1.4 


0.9 


1.9 


2.0 


1.4 


[N II] F 122 /F 205 


5.7 


3.7 


1.8 




6.2 


5.8 


4.1 


3.0 


3.7 


3.4 


2.5 


[N III] F57 Mm 


52 c 


55 


13 




115 


0.006 








62 


1.9 





[N III] F B7 /[N II] Fi 22 


21.0 


7.4 


5.0 




12.6 


0.0004 








9.1 


0.28 





[N III] F 57 /[0 III] F 52 


0.32 


0.60 


0.89 




0.33 


12 






0.53 


1.49 




[Ne III] F 36 /xm° 


17 








8.5 














[No II] F 12 8 ^m d 


1300 








1793 














[S IV] F10.5 Mm° 


100 








95 














[Ar III] F 9 .o A*m f 


350 








251 














[Ar III] F9.0 )im 6 


83 








150 














[Ar II] F7.0 )im s 


74 








13 














Ho + /H+ RRLs h 


0.045 








0.044 














[O III] F 52 ^m' 


194 








359 














[O III] Fss Mm' 


41 








56 














[N III] F57 Mm' 


50 








144 















b 



"Fluxes Xl0~ 18 W cm , corrected for extinction. 

b Model Abundances (xl0~ 6 ): C/H = 400, N/H = 125, O/H = 500, Nc/H = 340, S/H = 13.4, Ar/H = 4.5. 
c 0.75' Aperture, Colgan ct al. (1993). 
d 0.75' Aperture, Geballe et al. (1981). 
e 0.50' Aperture, Geballe ct al. (1981). 
f 0.67' Aperture, Geballe et al. (1981). 

g 0.35' Aperture, Cohen ct al. (1989), Simpson ct al. (1995b). 

h Radio recombination lines measured in 2.3' aperture, McGee & Newton (1981). 
'80" Aperture for ISO LWS. 
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Table 5 

W43 - Comparison of Models to Observations 



Observed Fluxes and Ratios a Model Fluxcs b and Ratios 



Line Flux or Ratio 


W43N 


W43C 


W43S 


Model D c 


Model S c 


Model T c 


Comments 




5.91 


2.25 


3.11 


3.50 


3.50 


3.50 




He+/H+ RRLs 




0.076 d 




0.075 


0.073 


0.075 


Good T eff 


[S III] Fig M m 


202 


122 


149 


155 


161 


163 




[S III] F33 Mm 


261 


149 


208 


149 


156 


160 




[S III] 19/33 fim 


0.78 


0.82 


0.72 


1.04 


1.03 


1.02 


Model JV e too high 


[O III] F52 Mm 


153 


130 


148 


205 


143 


81 




[O III] F ss pm 


59 


56 


64 


89 


62 


35 




[O III] 52/88 Mm 


2.6 


2.3 


2.3 


2.30 


2.33 


2.33 


Good N c 


[N II] F122 Mm° 


11 


4.1 


5.9 


4.2 


4.7 


4.2 




[N II] F 122 Mm 


4.9 


2.5 


2.6 


4.2 


4.7 


4.2 




[N II] F 205 Mm° 


1.14 


0.35 


0.42 


0.66 


0.75 


0.68 




[N II] 122/205 Mm 


4.3 


7.1 


6.2 


6.29 


6.28 


6.19 


Good N c 


[N III] F 57 Mm 


76 


61 


77 


61 


58 


65 




[N III] F 57 /[0 III] F 52 


0.50 


0.47 


0.52 


0.30 


0.41 


0.80 




[N III] F 57 /[N II] F122 


7 


15 


13 


15 


12 


15 




[N III] F 67 /[N II] F 122 


15 


23 


30 


15 


12 


15 




S++/H+ 


9.0 


14.1 


12.8 








Times 10~ 6 


S++/S 








0.919 


0.953 


0.965 


Implies S/H = 11.8, 11.4, or 11.3x10 


O++/S++ 


14.5 


19.5 


17.8 


0.583 f 


0.391 f 


0.221 f 


Implies O/H = 3.3, 4.8, or 8.3xl0~ 4 


N++/O++ 


0.521 


0.463 


0.516 


1.18 


1.58 


3.05 


Implies N/O = 0.42, 0.32, or 0.16 h 


N++/N+ 


0.6-2.5 


1.4-2.2 


1.3-3.2 


1.72 


1.45 


1.86 


Compare uniform Hc^/H^ 


N+/S++ 


3.0 


3.0 


2.6 


0.400 


0.428 


0.363 


Implies N/S = 7.25, 6.78, or 7.99 1 


N++/S++ 


7.5 


9.0 


9.2 


0.686 


0.620 


0.673 


Implies N/S = 12.55, 13.89, or 12. 79 1 



"Fluxes Xl0~ 18 W cm , corrected for extinction. Beam size is 45". 

b Model fluxes are scaled from the total model 5 GHz flux (~ 92 Jy) to 3.5 Jy since the KAO beams are much smaller than the total diameter 
of W43. 

c The model abundance ratios (xl0~ 6 ) are C/H = 400, N/H = 125, O/H = 500, Nc/H = 110, S/H = 12.1, and Ar/H = 4.5. 

d Weighted average of measurements of Churchwell et al. (1978), Lichten, Rodriguez, & Chaisson (1979), Lockman & Brown (1982), Peimbcrt 
et al. (1992). 

G Flux measured in a 60" beam and scaled to 45". 

f Model 0++/S++ is the density weighted ratio of <0++ /0>/<S+ + /S>, etc. (SCREH). 
g O/H = O/S x S/H. 

h N/H = 1.4, 1.5, or 1.4 xl0~ 4 if O/H = 3.3, 4.8, or 8.3 xl0~ 4 . 

'The average N/S = 9.9, 10.3, or 10.4, N/H = 1.17, 1.18, or 1.17 Xl0~ 4 , and O/H = 2.8, 3.7, or 7.2 Xl0~ 4 . 
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